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Abstract. We present the discovery of two Saturn-mass companions to HD 108147 and HD 168746. Both belong 
to the lightest ever discovered planets. The minimum mass of the companion to HD 168746 is of only 0.77 the 
mass of Saturn and its orbital period is 6.4 days. The companion to HD 108147 orbits its parent star in 10.9 days 
and its minimum mass is 1.34 that of Saturn. Its orbit is characterized by a high eccentricity, e = 0.50, indicating 
possibly the presence of a second companion. 

The detection of Saturn-mass planets by means of the Doppler technique demands high radial- velocity measure- 
ment precision. The two new candidates were discovered by means of the CORALIE echelle spectrograph. The 
instrumental accuracy of CORALIE combined with the simultaneous ThAr-reference technique has reached a 
level better than 3 ms _1 . On many observed objects the precision is now limited by photon noise. We present 
in this paper the weighted cross-correlation technique, which leads to an improvement in the photon noise of the 
computed radial velocity. We discuss as well a modification of the numerical cross-correlation mask which reduces 
significantly the residual perturbation effects produced by telluric absorption lines. 
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1. Introduction 

To date about 80 planetary candidates have been dis- 
covered. Their orbits show a widely-spread distribution 
of orbital parameters such as eccentricity and orbital pe- 



is large if its minimum mass rr^sini is large. This fact 
makes the technique particularly sensitive to high-mass 
companions. On the other hand, the detection of low-mass 
planets is more difficult and produces a detection bias on 



the low-mass end of the mass function (e.g. Qdry & Mayor 



ferent evolution mechanisms. The minimum mass of the 

detected planetary companions ranges from ~ 10 Jupiter 
masses down to about half the mass of Saturn. Four of 



riod (HjjB fc Mayor 2001p which might be related to dif- 20Q1 |. | Jorissenet aL gPOp, In order to reduce this bias the 



the known candidates have sub-saturnian masses ( Jorissen 



et al. 200 ip . The lightest planet, HD 83443 c, has been dis- 
covered by means of the CORALIE spectrograph and has 
a mass of only 0.53 M Sat flMayor et al~2000| ). 

The technique which underlies all these discoveries is 
that of precise radial-velocity measurements. It provides 
us with most orbital parameters. Because this technique 
does not allow us to constrain the angle of projection sin? 
of the orbital plane, only a minimum mass TO2 sini of the 
companion can be determined. The radial-velocity varia- 
tion in duced by a planetary companion on its parent star 



measurement precision must be improved by efficient and 
stable instrumentation and by efficient measurement and 
data reduction techniques. 

Using the CORALIE echelle spectrograph on the 1.2-m 
Euler Swiss telescope at LaSilla we are carrying out, since 
summer 1998, a large high-precision radial-velocity pro- 
gram ( iQueloz et al. 200C| ; |Udry et al. 2000p . Together with 
ELODIE in the northern hemisphere, CORALIE has al- 
lowed the discovery of about half of the known exoplanet 
candidates (Qdry & Mayor 2001). Recent improvement 



of the data reduction software have contributed to obtain 
an overall instrumental precision below ~ 3 ms _1 (Queloz 



et al. 2001a) over time scales longer than 3 years. The re- 
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* Based on observations collected with the CORALIE echelle 
spectrograph on the 1.2-m Euler Swiss telescope at LaSilla 
Observatory, ESO Chile 



& Carrier 2001) have proven that the short term precision 
is even better, namely about 1 ms -1 r.m.s over 1 night. 
The limitations are now determined mainly by photon 
noise and residual astroclimatic influence. 
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The present paper describes the discovery of two 
Saturn-mass companions to the stars HD 108147 and 
HD 168746. HD 168746 b, with its minimum mass of only 
0.77 Mg at , is one of the four sub-saturnian planets discov- 
ered to date. On the other hand, HD 108147 b distinguishes 
itself by a high eccentricity, fueling the discussion of the 
origin of such high eccentricity in the case of short-period 
companions. Before the orbital parameters and the stel- 
lar characteristics of these two objects described in detail 
in the second part of the paper, we discuss additional im- 
provements made recently in the data reduction, and more 
precisely, in the way of extracting the radial- velocity infor- 
mation obtained from the high-resolution spectra recorded 
with CORALIE. 

2. Improving the radial-velocity precision on 
CORALIE data 

2.1. The weighted cross-correlation 

The technique of the numerical cross correlation used to 
determine the precise radial velocity is described in detail 



in Baranne et al. (1996) and will not be discussed here. We 
only recall that the measured spectrum is correlated with a 
numerical mask consisting of 1 and value-zones, with the 
non-zero zones corresponding to the theoretical positions 
and widths of the stellar absorption lines at zero velocity. 
The cross-correlation function (CCF) is constructed by 
shifting the mask as a function of the Doppler velocity: 



CCF(vn) = J 5(A) • M(X VR )d\ (1) 
= J S(X) ■ ^ Mj(A„ R ) dX (2) 

i 

= J2 J S(\)M i (X VR )d\ = Y / CCF i (v K ), (3) 



where 




(4) 



In this equation S*(A) is the recorded spectrum, and 
M(X VR ) represents the Doppler-shiftcd numerical mask 
which can be expressed as the sum of masks Mi each 
corresponding to a stellar absorption line i. The result- 
ing CCF(vr) is a function describing somehow a flux- 
weighted "mean" profile of the stellar absorption lines 
transmitted by the mask. For the radial velocity value of 
the star, we take the minimum of the CCF, fitted with a 
Gaussian function. 

The cross-correlation technique has proven to be very 
robust and simple, and to deliver excellent results. The 
great advantage of the cross-correlation with a numeri- 
cal mask is that it does not need any high signal-to-noise 
reference spectra to compute the precise radial velocity. 



However, it has been shown in the past (Bouchy et al 



2001; Chclli 2000) that in terms of photon noise the tech- 



extract the radial-velocity information content in an op- 
timized way. Deep and sharp lines, for example, are not 
weighted sufficiently, although they contain more radial- 
velocity information than broad and weak lines. The chal- 
lenge was therefore to combine the simplicity and robust- 
ness of the cross correlation with a more efficient informa- 
tion extraction. 

The solution presented here assumes that the absorp- 
tion lines of the stellar spectrum used to compute the ra- 
dial velocity have all about the same FWHM, which is 
a good approximation since broad lines are not included 
in our standard numerical mask. The lines have different 
relative depths which are averaged in the resulting CCF. 
Lines with large relative depth contain intrisically more 
radial-velocity information than weak lines, however. Our 
goal is therefore to build up a weighted cross-correlation 
function CCF W which accounts for the correct weight of 
each spectral line contained in the mask: 



CCF™(v R ) 



CCF z (v R )- Wi 



(5) 



In order to understand the origin of uii the following ex- 
ample shall be presented: The contribution to the global 
CCF of a spectral line I with relative depth q = 0.5 
and a continuum level of Si = Iq is exactly the same 
as for a line m with relative depth c rn = 1 but a con- 
tinuum level at S m = Iq/2. The cross-correlation signal 
CCFi cx Ci ■ Si is indeed of equal amplitude in both cases, 
namely CCFi = CCF m cx Io/2. The noise on each point 
of the cross-correlation function is however proportional 
to the square root of the real spectral signal Si, and thus 
\/2 times larger in the CCFi. If we would fit each cross- 
correlation function CCFi independently we would get an 
precision on the fitted position which is V2 larger for the 
line m. Thus, the weight of the line I should be greater. 

In general, for a given amplitude of the cross- 
correlation signal CCFi, the noise on each point of the 
CCF - and thus on the resulting Gaussian fit - is of cx 
St cx CCFi cx — . The weight we have to give to the single 
CCFi must therefore be equal to u>.; = -ij = c^, where Ci is 

the relative depth of each spectral line i. Eq. (||) becomes 
then: 



^CCF 4 (« R )- Cl 

t 

J I(X)-J2(Mi(X v 



Ci) dX 



(6) 



(7) 



nique could still be improved. In fact, the CCF does not 



By means of a synthetic K0 dwarf spectrum we have com- 
puted the relative depth q (and thus the weight) for each 
absorption line contained in the numerical mask. We have 
modified the reduction software to take into account the 
relative weight of each line when computing the CCF. 
Finally, we have compared the old and the new code by 
numerical simulation on a set of 100 spectra issued from 
a high signal-to-noise spectrum on which statistical noise 
has been added. The result was a reduction of 1.25 of the 
measured radial velocity dispersion which is equivalent to 
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an hypothetical increase in signal-to-noise on the spec- 
trum of the same value. We have also applied the old and 
the new algorithm on the radial- velocity data of HD 108147 
and HD 168756. The results of this comparison will be pre- 
sented below. 



2.2. Telluric lines: reduced contamination effects 

Telluric absorption lines superimposed to the stellar spec- 
trum form absorption bands, especially in the red-visible 
wavelength region. The position of the telluric lines rela- 
tive to the stellar lines changes as a function of the relative 
radial velocity between Earth and star, while the strength 
of the telluric lines depends strongly on the astroclimatic 
conditions. When cross-correlating the recorded spectrum 
with the numerical mask a parasitic signal is produced by 
the casual correlation of a line in the mask with a telluric 
line in the spectrum. This might introduce an asymmetry 
in the CCF and consequently an error on the measured 
stellar radial velocity. In principle, the cross-correlation 
mask used in the ELODIE and CORALIE data reduction 
has been "cleaned" from telluric lines: all stellar lines for 
which the wavelength lies close (at a couple of line widths) 
to that of a telluric line have been removed from the nu- 
merical mask. Because of the quite large annual Earth mo- 
tion, however, the telluric lines change their relative posi- 
tion with regard to the stellar lines by about ±30kms _1 
(±0.61 at 6000 A). These means that every stellar line 
lying within ±30kms _1 from a telluric line is potentially 
affected by a telluric line more or less strongly in some 
period of the year. The superposition of these effects can 
then result in an increased "noise" or even in an erroneous 
periodic signal on the stellar velocity with a period of one 
year. Fig. [l] shows simulations made on real ELODIE spec- 
tra. 

In order to reduce this effect we have cleaned the origi- 
nal mask carefully from all stellar lines potentially affected 
by telluric lines. Although the photon-noise was slightly 
increased - since some regions of the stellar spectra were 
not used anymore in the cross correlation - quite impres- 
sive results have been obtained. In our simulations and as 
shown by the dashed curved in Fig. [l], the error produced 
by the telluric lines has been reduced at least by a factor 
of 3. 

Although in case of the CORALIE instrument the as- 
troclimatic conditions are in average more convenient than 
on ELODIE, and therefore the effect of the telluric lines 
is expected to be less important, we have decided to ap- 
ply the new numerical mask on the radial velocity data 
of HD 108147 and HD 168756 presented in this paper. The 
corresponding results are discussed below. 

Finally, we should also mention that the improvement 
of the mask presented in this section does not apply only 
to newly acquired data. It will allow us to re-correlate 
any spectrum obtained with ELODIE or CORALIE until 
the present and improve the existing radial-velocity data 
base. 
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Fig. 1. The plot shows the error on the measured radial 
velocity produced by telluric lines as a function of Earth 
radial velocity. For this purpose a real stellar spectrum has 
been Doppler-shifted by a radial-velocity value ranging 
from — 30kms _1 to -l-SOkms^ 1 and superimposed with a 
zero-velocity telluric absorption spectrum. The obtained 
radial velocity has then been compared to the star velocity 
obtained on the original stellar spectrum, the error being 
the difference between them. The solid line shows the re- 
sults using the original numerical mask while the dashed 
line represents the results with the new clean mask 

3. A Saturn-mass planet in eccentric orbit around 
HD 108147 

3.1. Stellar characteristics of HD 108147 

HD 108147 (HIP 60644) is a F8/G0 dwarf in the Crux 
constellation. Its magnitude is U = 6.99 while the 



HIPPARCOS catalogue ( |ESA 1997| ) lists a color index 
B — V — 0.537. The precise astrometric parallax is tt = 
25.95 ± 0.69 mas corresponding to a distance of about 
38.57 pc from the Sun. The derived absolute magnitude, 
My = 4.06, is typical for a GO dwarf. 

Stellar parameters such as effective temperature 
T e g = 6265 K, surface gravity logg = 4.59, as well as 
[Fe/H] = +0.2 have been derived in the detailed LTE spec- 
troscopic analysis carried out by Santos et al. (2001). The 
obtained metallicity is slightly higher than the average 
value for stars of the CORALIE sample, like most of the 
stars with giant planets, and is very close to the mean 



value of [Fc/H] of stars with planets (Santos et al. 2001). 
Using the evolutionary tracks of the Geneva models given 



by Schaerer et al. (1993), Santos et al. (2001) compute a 
stellar mass M = 1.27M Q . The stellar mass is higher than 
the "typical" mass of GO dwarfs and can be explained by 
the high metallicity of the star. 

Fig. U shows the CaiiH absorption line region of the 
CORALIE spectrum at A3968.5A. The emission flux in 
the core of the Ca n H line corrected for the photospheric 
flux provides us with the chromospheric activity index 
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Ca II H 




Table 1. Observed and inferred stellar parameters for 
HD 108147 and HD 168746. Photometric and astrometric 
data have been extracted from the HIPPARCOS catalogue 



3964 



3966 



3968 3970 

\(A) 



3972 
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Fig. 2. A3968.5A CaiiH absorption line region of the 
summed CORALIE spectra. Upper panel: HD 108147. 
Lower Panel: HD 168746. Both objects show low chromo- 
spheric activity. Strong activity would result in an emis- 
sion peak in the center of the absorption line 



Scor ( [Santos et al. 200~ ) from which we derive the activ- 
ity indicator log(R' liK ) = —4.72. This value is typical for 



stars with a low chromospheric activity level ( Henry ct al 
1996()7]U sing the calibration given by ponahue (1993|) and 



quoted in |Hcnry ct al. (1996| ) we compute for this star an 
age of approximately 2 Gyrs, while th e rotational period 
resulting from the calibration given by Noyes et al. (1984 ) 
is of 8.7 days. The star is not seen as photometrically vari- 
able in the HIPPARCOS data, confirming again the low 
activity level. We have estimated the projected rotational 
velocity of the star to be i»sini = 5.3 kms" 1 by means 
of the CORA LIE cross-correlation function (CCF) ( |Queloz 



et al. 1998 ). The relatively high stellar rotation could 
cause a small jitter on the radial-velocity data, which are 
however expected to be in the order of only few ms" 1 
( fcaar fc Donahue 1997| ; |Santos et al. 200 1| ). The observed 
and inferred stellar parameters are summarized in Table |l|. 

3.2. CORALIE orbital solution for HD 108147 

The precise radial velocity data of HD 108147 have 
been collected by our group during the period of time 
from March 1999 to February 2002[|. The result of this 
campaign is a set of 118 data points having a mean 
photon-noise error on the individual measurements of 
(si) =7. 8ms -1 . A periodic variation of the radial velocity 
could be detected on this data, which clearly indicates the 



(ESA 1997) while spectroscopic data are from Santos et al 
(200 ID 



Parameter 


HD 108147 


HD 168746 


Spectral Type 


F8/G0 


G5 


V 


6.99 


7.95 


B — V 


0.537 


0.713 


7r |masj 


25.95 ± 0.69 


23.19 ± 0.96 


Distance [pc] 


OO _1_ 1 

38.0/ ± 1 


4o.l2 ± 1.0 


M v 


4.06 


4.78 


L/L Q 


1.93 


1.10 


[Fe/H] 


0.2 ±0.06 


-0.06 ±0.05 


M/Mq 


1.27 ±0.02 


0.88 ±0.01 


Toff [K] 


6265 ± 40 


5610 ± 30 


log g [cgs] 


4.59 ±0.15 


4.50 ±0.15 


vs'mi [kms -1 ] 


5.3 


1.0 


l0g(^ HK ) 


-4.72 




P ro t(i*H K )[days] 


8.7 




age(i& K ) [Gyr] 


2.17 


-/old star 



presence of a planetary companion. The measured varia- 
tion cannot be produced by stellar activity: The Geneva 
photometry data show very low dispersion of 1 m mag and 
no CCF-bisector variation ( Qucloz ct al. 2001b ) has been 
measured at the ms -1 level. Therefore the planetary ex- 
planation seems to be the most likely. 

The best-fit Keplerian orbit to the data is shown 
in Fig. It yields a precisely-determined orbital pe- 
riod P of 10.901 ± 0.001 days and a large eccentricity 
e = 0.50 ± 0.03. The semi-amplitude of the radial- velocity 
variation is K = 36 ± 1 ms -1 . The weighted r.m.s. of the 
data to the Keplerian fit is 9.2 ms" 1 . The complete set 
of orbital elements with their uncertainties are given in 
Table |. 

Using the best-fit orbital parameters and the mass of 
HD 108147 given above we derive for the companion a min- 
imum mass ui2 smi = 0.40 Mj up (which is about only 
1.34 times the mass of Saturn). Because of the many data 
points and the long observation period the orbit is deter- 
mined very accurately. This allows us to determine the 
minimum mass of the companion with accuracy of bet- 
ter than 4%, provided that we do not consider the major 
error source, namely the uncertainty on the mass of the 
primary. From the orbital parameters and the star mass 
we get also the separation of the companion to its parent 
star which is a = 0.104 AU. The surface equilibrium tem- 
perature of the planet at such a distance is estimated to 



1 The discovery was announced in May 2000, 
www.eso.org/ outreach / press-rel / pr-2000 / pr-13-00.html 



be about 890 K, following |Guillot et al. (L996| ). 

HD 108147b belongs to the so-called hot Jupiter (or 
better: hot Saturn!) category of extra-solar planets. The 
close location to its parent star makes the planet a good 
candidate for a photometric transit search. The photomet- 
ric monitoring described in a forthcoming paper (Olsen et 
al., in prep.) did unfortunately not show any indication 
for a transit. 
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HD 108147 




JD-2400000 [days] 

Fig. 3. Phase-folded radial-velocity measurements ob- 
tained with COR AT, IE for HD 1 081 47. The error bars rep- 



resent Iphoton-noisc errors only. On the lower panel the 
residuals of the measured radial velocities to the fitted or- 
bit are plotted as a function of time. They show a tiny 
indication for the presence of a long-period, second com- 
panion of HD 108147 

Table 2. CORALIE best Keplerian orbital solutions de- 
rived for HD 108147 and HD 168746, as well as inferred 
planetary parameters 



Parameter 


HD 108147 


HD 168746 


P 


[days] 


10.901 ±0.001 


6.403 ±0.001 


T 


[JD] 


2451591.6 ±0.1 


2451994.7 ±0.4 


e 




0.498 ±0.025 


0.081 ±0.029 


7 


[kms" 1 ] 


5.026 ±0.001 


-25.562 ±0.001 


Ul 


[deg] 


318.95 ±3.03 


16.30 ±20.88 


K 


[ms- 1 ] 


36±1 


27±1 


-/Vmc as 




118 


154 


a(O-C) 


[ms- 1 ] 


9.2 


9.8 


ai sin i 


[Mm] 


4.634 


2.364 


f(m) 


[M ] 


3.337-10" 9 


1.284-10" 9 


mi 


[M ] 


1.27 


0.88 


iri2 sini 


[Mj U p] 


0.40 


0.23 


a 


[AU] 


0.104 


0.065 


T c q 


[K] 


890 


900 



3.3. A planetary system around HD 108147 ? 

The weighted r.m.s. of the data to the Keplerian fit is 
9.2 ms" 1 and the reduced \ 1S 1-440, indicating an inter- 
nal error for the single measurements of 6.4 ms -1 . The 
residual dispersion is 6.6 ms" 1 and is only partially com- 
posed of instrumental errors, which are in the order of 
3 ms" 1 (Qucloz ct al. 2001a). The remaining 6ms" 1 could 



Table 3. Comparison of the different data reductions on 
HD 108147. The r.m.s of the data to the Keplerian fit is 
reduced from 11.9 to 9.2 ms" 1 if the weighted cross cor- 
relation and the clean numerical mask are used 



Cross correlation 


Obtained r.m.s [ms A ] 


Standard 


11.9 


Clean mask only 


10.7 


Weighted only 


10.1 


Weighted & clean mask 


9.2 



fating. Closer analysis of the residuals (Fig. [§[ bottom) 
might indicate however the presence of a possible second 
companion with an orbital period P of 587 ± 34 days and 
a minimum mass of mi sini = 0.42 Mj up . A best-fit solu- 
tion considering a second companion reduces the r.m.s to 
the Keplerian to 8.0 m s" 1 while the reduced x passes from 
1.440 to 1.290. In order to confirm or reject this possibility 
we will continue to collect additional high-precision data 
during the next observational seasons. We plan to moni- 
tor the radial velocity of this object using the CORALIE 
spectrograph, but also taking advantage of the superior 
performances of the HARPS^ spectrograph (Pepe ct al 



2000) which is expected to reach a radial- velocity preci- 



sion of 1 m s . 

3.4. Results of the improved cross correlation 

As mentioned previously, the radial velocity data have 
been computed by using the weighted cross correlation 
and the new numerical mask. For comparison we have 
nevertheless computed the radial velocity data also us- 
ing the standard on-line data reduction. The results are 
summarized in Table |^. 

In the error budget of our best single-orbit solution 
about 6.6ms" 1 arise from sources other that photon 
noise. The r.m.s obtained with the standard algorithm 
was 11.9ms" 1 and with the weighted cross correlation 
10.1 ms" 1 . Considered this, the weighted cross correlation 
reduced the radial velocity dispersion arising from photon 
noise by a factor 1.31, confirming the simulation results. 
On the other hand, the r.m.s of 11.9 ms" 1 obtained with 
the original mask passes to 10.7ms" 1 with the new clean 
mask. Thus, we estimate that the telluric lines add a dis- 
persion on the stellar radial velocity of roughly 5ms" 1 . 

3.5. The eccentricity of the orbit around HD 108147 

Considered the short orbital period, the eccentricity of the 
orbit around HD 108147 is surprisingly high. According to 
the commonly believed migration scenario for the forma- 
tion of hot Jupiters, the giant planets like HD 108147 can- 
not be formed in situ at a distance of only 0.1 AU from 
the parent star. On the other hand the orbit of planets 
having suffered a strong migration are expected to be cir- 



be due to stellar jitter, since the star is fairly rapidly ro- 



2 HARPS will be installed on the ESO 3.6-m telescope and 
become operational at the beginning of 2003 
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cularized by the tidal interaction with the accretion disk 
flGoldrcich fcTV cmainc 1980| ). Another mechanism must 
therefore be at the origin of its high eccentricity. One suit- 
able explanation is provided by the gravitational interac- 
tion between multiple giant planets in the early stages 
of the system formation ( Weidenschilling & Marzari 1996; 
Rasio & Ford 1996; Lin & Ida 1997). As a result one planet 



HD 168746 



Coralie 



could have been projected on a short-period orbit with a 
high eccentricity, leaving different possibilities concerning 
the destiny of its scattering partner. Eccentric orbits can 
however also be produced by additional companions on 
longer-period orbits, which perturb the short-period or- 



bit of the inner companion by tidal forces ( Mazeh et al 
19970The companion can be either planetary or stellar 
and would be seen in many cases only as a slow drift of the 
gamma-velocity 7. Finally, additional scenarios to explain 
high eccentricities were proposed recently: For example 
the migration of a single Jupiter-mass planet through a 
planetesimal disk ( Lcvison et al. 1998 ), or, the simulta- 
neous formation and migration of two or more planets 
through the accretion disk, e.g locked in a resonant sys- 
tem ( [Murray et al. 2002| ). The numerical simulations in 
the latter paper show that a large variety of possible re- 
sults can be obtained, depending on the initial planetary 
mass(es), disk properties and time scales. 

In summary, various mechanisms can explain short- 
period eccentric orbits. In almost all of them interaction 
with planetary or stellar companions plays a fundamental 
role, and the large eccentricity of the orbit might be an 
indication for their presence. Therefore it is very likely to 
find a second, long-period companion around HD 108147. 
As mentioned above the long-term follow up of this object 
might confirm or reject this possibility. 

4. A Sub-Saturn mass planet around HD 168746 

4.1. Stellar characteristics of HD 168746 



HD 168746 (HIP 90004) is a G5 dwarf of magnitude 
V = 7.95 located at the boundary between the Scutum and 
the Serpens Cauda constellations. The HIPPARCOS cat- 
alogue lists a color index B — V = 0.713. The catalogue 
indicates also an astrometric parallax of tt = 23. 19 mas 
( ESA 1997 ) corresponding to a distance of about 38.57 pc 
between the star and the Sun. The resulting absolute mag- 
nitude is M v = 4.78. 



The spectroscopic analysis by Santos et al. (2001) indi- 
cates an effective temperature T e g = 5610 K, log<? = 4.50, 
and [Fe/H] = —0.06. The combination of the precise spec- 
troscopic parameters with the evolution models provide a 
stellar mass of M = 0.88 M@. The emission flux in the 
core of the Ca 11 H shows no evidence for chromospheric 
activity (see Fig. |2|). Since the star is too faint it has not 
been possible to compute a reliable value for i?Jj K from the 
CORALIE spectra. Therefore it was not possible to derive 
neither the age nor P ro t- Evolutionary tracks suggest for 
this star an age of at least several Gyrs. 



^25. 55 




51400 51600 51800 52000 52200 

JD-2400000 [days] 

Fig. 4. Phase-folded radial-velocity measurements ob- 
tained with CORALIE for HD 168746. The error bars rep- 
resent photon-noise errors only. The lower panel shows the 
residuals of the measured radial velocities to the fitted or- 
bit as a function of time 



4.2. CORALIE orbital solution for HD 168746 

During the period spanning from May 1999 to September 
2001^ we have recorded 154 precise radial-velocity data 
points of HD 168746. Because of the magnitude of the star 
the mean photon-noise error on the individual measure- 
ments is relatively high, namely of (£j) =9.7ms . Due 
to the large amount of data we could nevertheless iden- 
tify easily the presence of a companion with a very low 
minimum mass. Because of the photometric stability of 
HD 168746 and the absence of line-bisector variations we 
can exclude that the measured radial- velocity variation is 
produced by stellar activity. 

Fig. |] shows the best fitted Keplerian orbit for the 
companion of HD 168746. We deduce an orbital period P 
of 6.403 ± 0.001 days and a rather small eccentricity of 
e = 0.08 ±0.03. The semi-amplitude of the radial-velocity 
variation is only K — 27 ± lms . The period and the 
semi-amplitude are well constrained thanks to the large 
amount of data points covering many periods. The com- 
plete set of orbital elements with their uncertainties are 
given in Table 

The best-fit orbital parameters and the mass of 
HD 168746 yield for its companion a minimum mass of 
mi sini = 0.23 Mj up , which is only 0.77 times the mass 
of Saturn. The accuracy obtained on this mass estimation 
is better than 3%, if we again do not consider the uncer- 
tainty on the mass of the primary. The separation of the 



The discovery was announced in May 2000, 
www.eso.org/outreach/press-rel/pr-2000/pr-13-00.html 
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companion to its parent star is a = 0.065 AU. The surface 
equilibrium temperature of the planet at such a distance 
is about 900 K. Similarly to HD 108147, a dedicated search 
for a possible photometric transit of the planets in front of 
HD 168746 remained unsuccessful (Olsen et. al, in prep.). 

The radial velocity data were computed, as for 
HD 108147, using the weighted correlation and the cleaned 
numerical mask. We confirm a reduction of the mea- 
sured dispersion: The weighted r.m.s. of the data to the 
Keplerian fit decreased from the original 11.2ms _1 to the 
final 9.8ms" 1 , with a reduced \ °f 1-5. 

The CORALIE individual radial-velocity measure- 
ments presented in this paper are available in electronic 
form at CDS in Strasbourg. 

5. Concluding remarks 

The companions to HD 108147 and HD 168746 possess 
minimum masses of 1.34 and 0.77Mg at , respectively. The 
latter is one of the lightest extra-solar planet candidates 
found to date. Both objects have been detected by means 
of the CORALIE spectrograph. CORALIE and its twin 
instrument ELODIE are particularly efficient for the de- 
tection of light planets: Seven of the ten lightest candi- 
dates and two of the four sub-saturnian companions have 
been discovered using these instruments^, among which 
HD 83443 c, the lightest ever discovered planetary com- 
panion with only 0.53 Ms at . 

The suitability of CORALIE for the detection of low- 
mass planetary companions arises from the intrinsic in- 
strument precision which is better than 3ms _1 over long 
time scales. Despite the modest size of the telescope the 
overall detection efficiency of CORALIE is high, due cer- 
tainly to the telescope-time availability and, in particular, 
to the simultaneous ThAr reference. The ThAr technique 
allows us to cover the full visible wavelength range with- 
out any transmission losses, ensuring a gain in flux of at 
least a factor of 6 compared to the use of an iodine absorp- 
tion cell. The photon noise is thus reduced by a factor of 
about 2-2.5 compared to an exposure of same duration us- 
ing the iodine cell (Bouchy et al. 2001). Nevertheless, the 
measurements remain in many cases photon-noise limited. 
The new weighted cross-correlation algorithm presented in 
this paper allowed us to reduce by about 1.25 the photon- 
noise contribution to the radial-velocity data. Although 
some improvement of CORALIE and its data reduction 
software is still going on, our major future step will be to 
complete the realization of HARPS (Pepe et al. 200C), a 
new stable and efficient echelle spectrograph dedicated to 
the search for extra-solar planets with the ESO 3.6-m tele- 
scope at La Silla. The aimed precision of 1ms -1 should 
considerably increase the detection rate of sub-saturnian 
planets. 
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